Metallic "Ferroelectricity" in the Pyrochlore Cd2Re207 
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A class of materials known as "ferroelectric metals" was discussed theoretically by Anderson and 
Blount in 1965 [Phys. Rev. Lett. 14, 217 (1965)], but to date no examples of this class have been 
reported. Here we present measurements of the elastic moduli of Cd2Re2 07 through the 200 K 
cubic-to-tetragonal phase transition. A Landau analysis of the moduli reveals that the transition is 
consistent with Cd2Re207 being classified as a "ferroelectric metal" in the weaker sense described 
by Anderson and Blount (loss of a center of symmetry). First-principles calculations of the lattice 
instabilities indicate that the dominant lattice instability corresponds to a two-fold degenerate mode 
with Eu symmetry, and that motions of the O ions forming the O octahedra dominate the energetics 
of the transition. 
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Rinerant electrons screen electric fields and inhibit the 
electrostatic forces responsible for ferroelectric distor- 
tions. Therefore, in a metallic system one does not ex- 
pect to find structural transitions similar to those found 
in insulating materials with a tendency toward ferroelec- 
tricity. The idea that metallic behavior and ferroelectric- 
ity may not always be incompatible had an early cham- 
pion in B. T. Matthias 0. A groundbreaking paper on 
this subject was written by P. W. Anderson and E. I. 
Blount (A&B) in 1965 2]. Applying Landau theory to 
a continuous cubic-to-tetragonal (C-T) structural phase 
transition, A&B concluded that "a transition from cu- 
bic to tetragonal in which the only order parameter in 
Landau's sense is the unit cell shape, i.e., the strain, 
can be second order only with probability zero." Apply- 
ing these ideas to the C-T transition found in A- 15 su- 
perconductors such as VaSi and NbsSn, A&B concluded 
that "these and perhaps several other metallic transitions 
may be 'ferroelectric' in the sense of the appearance of 
a polar axis, or possibly at least involve the loss of a 
center of symmetry." Ultimately, however, it was shown 
that the structural transitions in the A- 15 compounds 
were not continuous but rather weakly first order, and 
that strain was indeed the appropriate order parameter 
0. As no other materials seemed to fit A&B's criteria, 
ideas regarding metallic "ferroelectricity" have not been 
pursued over the past thirty-eight years. 

The pyrochlore Cd2Re207 has attracted attention re- 
cently as an oxide superconductor on a geometrically 
frustrated lattice 0, IE IE S 13 • The normal state prop- 
erties of Cd2Re2 07 are also intriguing [E 0, O, 0, 
EE Q m, particularly the C-T transition at 200 K 
that profoundly affects the electronic structure, trans- 



port, and magnetic susceptibility of this material. Ev- 
idence from resistivity, specific heat, NQR, and X-ray 
diffraction |E EE EE EE 113 indicates that the 200 K 
C-T transition is continuous. Another, first-order, struc- 
tural phase transition (Tc = 120 K ) has been reported 
in Cd2Re207 EEl; but because the issues raised by this 
transition are secondary to the main focus of this pa- 
per, we will only mention this lower-temperature tran- 
sition briefly in what follows. Despite several studies 
EE EE El, the structure of Cd2Re207 below 200 K has 
not been fully determined, mainly because the depar- 
ture from cubic symmetry is extremely small and even 
high-resolution X-ray measurements [la, |l3| can barely 
detect the splitting of the cubic Bragg peaks. It has 
been shown, however, that there is no multiplication of 
the unit cell below the 200 K transition (the transition 
is ferrodistortive) , and that a loss of three- fold symme- 
try accompanies the transition 0- There is also some 
evidence that inversion symmetry may be broken in the 
low-temperature phases [UI. Cd2Re207, therefore, is a 
good candidate for becoming the first material to obey 
A&B's criteria provided that strain can be ruled out as 
the primary order parameter for the C-T transition at 
200 K. 

Single crystals of Cd2Re207 were grown from the va- 
por in sealed silica tubes using Cd metal (5]X^^ Johnson 
Matthey) and Re207 (3N, Johnson Matthey) [EIE|. The 
measurements reported in this paper were performed on a 
crystal cut into a rectangular parallelepiped with dimen- 
sions 1. 2x1. 9x2. 3 mm'^. The sample was oriented with 
all faces perpendicular to the crystallographic < 100 > 
axes. The experimental density of the sample was 8.795 
g/cm?; this can be compared with the x-ray density of 
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FIG. 1: Elastic moduli vs. temperature obtained on a single 
crystal of Cd2Re207 using RUS. The relationship between 
sound velocity and elastic moduli is given by Va = (C/p)^^'^, 
where C is the eflFective elastic constant and p is the density. 
In a cubic crystal, longitudinal waves propagating in the [100] 
direction are governed by Cn, and both transverse waves are 
governed by C44. In the [110] direction, longitudinal waves are 
governed by Ci[iio] = l/2(Cii + C12 + 2C44), one transverse 
wave is governed by C44, and the other by Ct[iioi ~ l/2(Cii — 
C12). In the [111] direction, longitudinal waves are governed 
by C^fiii] = l/3(Cii + 2Ci2 + 4C44), and both transverse 
waves are governed by Ct[iii] = l/3(Cii — C12 + C44). Note 
the small magnitude (2%) of the anomaly in C44 compared 
to the anomalies in the other moduli. 



8.814 g/cm^. 

Resonant Ultrasound Spectroscopy (RUS) measure- 



ments were performed as a function of temperature (5- 
300 K) to determine the elastic moduli of the sample. 
RUS is a technique developed by Migliori, et al. for 
determining the complete elastic tensor of a small sin- 
gle crystal by measuring its free-body resonances. This 
method has the advantage that all moduli can be de- 
termined simultaneously, thereby avoiding remounts of 
transducers and multiple temperature sweeps. 

In Fig. 1 we plot the elastic moduh of Cd2Re207 vs. 
temperature for the three modes of elastic waves in the 
principal propagation directions in a cubic system [20l [. 
Deep into the transition the ultrasonic absorption of the 
sample became so great that for several temperatures not 
enough resonances were observed to allow for an accurate 
determination of all three elastic moduli. However, the 
lowest frequency resonance depended almost exclusively 
on C44 and was visible throughout the entire transition. 

If strain were the order parameter, we would expect 
at least one of these elastic constants to soften dramat- 
ically as expected for an elastic instability. This is not 
observed. The salient feature of Fig. 1 is the step-like 
change at 200 K in all of the moduli except C44. This 
behavior of the moduli can be modeled using the order 
parameter proposed by Sergienko and Curnoe This 
order parameter involves collective atomic displacements 
corresponding to a long-wavelength phonon of sym- 
metry. 

The minimal model for Landau free energy F which 
accounts for the anomalies of the elastic moduli should 
include a ferrodistortive energy expanded in terms of 
the structural order parameter (ryi, 7^2), the elastic energy 
Fei , and coupling between the order parameter and strain 
Fd-pi:- 



F,i = l/2C0i(e2 + e2 + ^2) ^ ^02(6162 + 6263 + eieg) + 1/20^4 + + ei), 
Fd-ei = Xiivl + V^){ei + 62 + 63) + X2[{vf - V^){ei + 62 - 263) + 2^/8771772(61 - 62)] 

+^H{vl + + 4 + 4) + ^^2[{r^l - vDi^l + 2ei) + 2^3771772(6! - ei)]. 



with F = Fd + F^i+ Fd^ei- Here Cfi, C^a, and CJ4 are 
the elastic moduli in the cubic phase (771 = 772 = 0) and 
01,02,03, and 61 are the Landau expansion coefhcients, 
ei is the strain, and Xi and /i^ are coupling constants. 

Terms up to fourth order in Fd are isotropic, there- 
fore sixth order terms are required to lift the degener- 
acy between ordered states. Since the phase transition 
at Tc — 200 K is second order, the coefficient 02 must 
be positive. The elastic moduli in the tetragonal phase 
(771 = 0, 7/2 7^ 0) can be calculated using Slonczewski- 



Thomas formalism 22] 

^ d^F d^F d^F /d^Fy^ 

""^ dcidej 9772 9ei dri2dej \ dr^l J 

where the equibrium values of 772 and the strain tensor 
are calculated from the system of equations 

1^ = 1^ = 0. . = !,.. ..6. (2) 

0772 oei 
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where ^ = Ai — A2, B = Ai4-2A2. Higher order terms can 
easily be calculated from the above equations by expand- 
ing in 772 but the resulting expressions are cumbersome. 
Steps are therefore expected in the elastic moduli Cn, 
C33, C12, and C13 at the continuous transition, while 
the shear moduli C44 and Cge have continuous anoma- 
lies. The data in Fig . 1 are qualitatively consistent with 
these predictions |23l |. 

In many respects the elastic behavior of Cd2Re2 07 
(CRO) resembles that of SrTiOa (STO), but there are 
some important differences. In STO the coupling is lin- 
ear in all components of the strain but quadratic in the 
order parameter 0, IH> 1^ Ell , whereas in CRO it is 
linear only in the diagonal components ei, 62, 63 of the 
strain tensor. It should also be kept in mind that STO is 
antiferrodistortive and inversion symmetry is not broken 
in the tetragonal phase. Like in STO, the measured elas- 
tic anomalies of Cn and (Cn — Ci2)/2 in CRO are not 
true step functions but are broadened by several degrees. 
Liithi and Moran ascribe this behavior in STO to 
residual strain, and there is no reason not to expect this 
strain in CRO as well. The downward "dip" observed 
near Tq in the longitudinal elastic moduli has also been 
observed in STO, and can be ascribed to order param- 
eter fluctuations |2.'tI |. Domain formation below Tc in- 
troduces several difficulties into a quantitative analysis 
of the data, both because of the anisotropics associated 
with domains and because domain wall motions make 
important contributions to the elastic moduli. To obtain 
elastic moduli below the transition we assumed a ran- 
dom domain distribution and that the sample retained a 
macroscopic cubic symmetry. Resolving these difficulties 
requires a detailed understanding of the microstructure of 
Cd2Re207, and this knowledge is not yet available. How- 
ever, none of these difficulties affects the phenomenology 
presented here. 

To help elucidate the origin of the lattice instabilities, 
first principles calculations were performed in the local 
density approximation (LDA) using the general potential 
linearized augmented planewave (LAPW) method jl^ 
as described in Ref . |27l l28l | . Initially, scalar relativis- 
tic calculations of the atomic forces were performed for 
a sufficient number of small atomic displacements away 
from the equilibrium structure to determine the dynami- 
cal matrix for an 88 atom supercell. Then the full phonon 
dispersion relations were obtained by a direct method 
using the PHONON program [2^. At the zone center 
two very unstable modes were found, a two-fold degen- 
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FIG. 2: Relativistic LDA energetics of lattice instabilities in 
Cd2Re207. The coordinates of the 0(1) atoms are in Carte- 
sian coordinates, units of the lattice constant a = 10.219 A: 
(0.315- 5e + 5t, 0.625, 0.625), (0.935 + 5e + St, 0.625, 
0.625), (0.625, 0.315, 0.625), (0.625, 0.935, 0.625), (0.625, 
0.625, 0.315 + Se), (0.625, 0.625, 0.935 - Se), (0.375, 0.065, 
0.375), (0.375, 0.685, 0.375), (0.065 + Se + St, 0.375, 0.375), 
(0.685 - Se + St, 0.375, 0.375), (0.375, 0.375, 0.685 + Se), 
(0.375, 0.375, 0.065 - Se), where E and T are the Eu and Tin 
frozen-in amplitudes, respectively. 



erate Eu symmetry mode, and a three-fold degenerate 
Tiu mode. Both of these modes involved breaking of in- 
version symmetry and have eigenvectors that are heavily 
dominated by the 0(1) site (forming the O octahedra co- 
ordinating the Re ions). In addition, two more weakly 
unstable and several low frequency but stable modes in- 
volving Cd and Re atom shifts were found, many of them 
also breaking inversion symmetry. 

Since the unstable modes were heavily dominated by 
0(1) motion, we neglected the minor components and 
calculated the energetics as a function of distortion am- 
plitude with a tetragonal cell for the Eu and Tiu displace- 
ment patterns. These calculations were done relativisti- 
cally, including spin-orbit, which we find to significantly 
affect the energetics by reducing the tendency towards 
lattice instability. Nonetheless, we still find a substantial 
instability of the Eu mode and a marginal instability of 
the Tiu mode, as shown in Fig. 2. Clearly, the dominant 
instability corresponds to the Eu mode, and so we expect 
that the higher temperature phase transition is due to its 
freezing in. 

The lower temperature transition might be associ- 
ated with a change of symmetry into which this mode 
freezes (note that it is two-fold degenerate), as in, e.g., 
the cubic-tetragonal-orthorhombic-rhombohedral transi- 
tion sequence of BaTiOa under cooling. Freezing of the 
doubly degenerate Eu mode alone may result in three 
possible low symmetry structures /4m2, /4i22, and F222 
|21|. We have performed the LDA calculations for these 
three space groups and found that the energetics are 
nearly the same. This can be understood because terms 
up to fourth-order in are isotropic. In Fig. 2 we show 
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the calculations, corresponding to the /4i22 structure, 
which was proposed to be the lowest temperature phase 
in Ref. 0|. In this case, soft modes associated with the 
could play an important role in the low-temperature 
superconductivity. 

Alternately, we note that the two unstable modes are 
of different symmetry, and therefore do not interact at 
lowest order. Perhaps a mode related to the Ti„ 0(1) 
displacement but with additional metal and/or 0(2) dis- 
placement is unstable enough to freeze in and give the 
lower temperature transition. 

The temperature dependence of the elastic moduli 
shown in Fig. 1 allows us to rule out strain as an order 
parameter in Cd2Re207. The conclusions of Anderson 
and Blount can therefore be applied to Cd2Re2 07: 
the order parameter must either be (1) "some electronic 
mystery parameter," or (2) "some change in symmetry, 
such as the loss of the inversion center." At present, the 
evidence supports the second possibility with most likely 
candidate for the order parameter being a small, coher- 
ent (hence "ferroelectric" ) collective displacement of the 
atoms with symmetry, dominated by the motion of 
the 0(1) atom as discussed above. Even though, fol- 
lowing A&B, we use the term "ferroelectric," we would 
like to stress that although inversion symmetry is broken 
there is no evidence that a polar axis is formed; indeed, 
our analysis is full y co nsistent with the /4m2 space group 
proposed in Ref. |l7l |. and this space group is piezoelec- 
tric, although because Cd2Re207 is metallic, no piezo- 
electric coupling is possible here. 

It is an open question whether the concept of a 
"ferroelectric metal" will be fruitful in explaining the 
anomalous physical properties at the 200 K transition 
in Cd2Re207. In many ways, the 200 K transition is 
reminiscent of a charge density wave (CDW) transition 
such as that observed in TaSe2, although CDWs in cu- 
bic materials such as pyrochlores are generally not ex- 
pected. Moreover, electronic structure calculations in- 
dicate a nearly isotropic Fermi surface and no obvious 
nesting or CDW instability 26] . If Cd2Re207 is indeed 
a "ferroelectric metal" however, one can imagine a re- 
distribution of charge within the material and physical 
properties that mimic a CDW transition. Also, the dra- 
matic decrease of the electrical resistivity of Cd2Re207 
below 200 K finds a natural explanation in terms of re- 
duced scattering from the unstable ions as they freeze in. 
It is hoped that the identification of Cd2Re2 07 as a "fer- 
roelectric metal" will stimulate theoretical development 
of the unique continuous phase transition found in this 
material. 
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